Protein and antigen profiles of third-stage larvae of Gnathostoma spinigerum assessed with next-generation sequencing transcriptomic information

Gnathostomiasis is a food-borne zoonotic disease that can affect humans who eat improperly cooked meat containg infective third-stage larvae. Definitive diagnosis is through larval recovery. However, this is an invasive technique and is impractical if the larvae have encysted in inaccessible areas of the body. Antigen or antibody detection might be more interesting techniques for diagnosis. Proteomic could elucidate diagnostic markers and improve our understanding of parasite biology. However, proteomic studies on Gnathostoma spinigerum are hampered by the lack of a comprehensive database for protein identification. This study aimed to explore the protein and antigen profiles of advanced third-stage G. spinigerum larvae (aL3Gs) using interrogation of mass spectrometry data and an in-house transcriptomic database for protein identification. Immunoproteomic analysis found 74 proteins in 24-kDa SDS-PAGE bands, which is size-specific for the immunodiagnosis of gnathostomiasis. Moreover, 13 proteins were found in 2-DE 24-kDa bands. The data suggest that collagenase 3, cathepsin B, glutathione S-transferase 1, cuticle collagen 14, major antigen, zinc metalloproteinase nas-4, major egg antigen, peroxiredoxin, and superoxide dismutase [Cu–Zn] may be good candidates for novel human gnathostomiasis diagnostic assays. These findings improve our understanding of the parasite’s biology and provide additional potential targets for novel therapeutics, diagnostics, and vaccines.

Typically, aL3Gs are harvested from eel livers. The process of CWA production is complex and complicated by the seasonal prevalence of the aL3Gs in eels 7,8 . With advances in proteomics and mass spectrometry, the 24-kDa antigen has been further characterized and found to have high homology to peptide sequence regions of cyclophilin, actin, matrix metalloproteinase-like protein, and intermediate filament protein B. Related studies have shown the applicability of some of these peptides as antigens 9,10 . However, no commercial diagnostic tests are currently available for gnathostomiasis.
For many years, there was no effective treatment for gnathostomiasis and currently the only successful treatment option is surgical excision of the larvae. Various medications, including thiabendazole, praziquantel, metronidazole, diethylcarbamazine, and quinine, have been explored in animal models and humans without success 11 . Proteomic technology has provided crucial information about cellular and molecular processes, leading to an improved understanding of the biology of many different parasites. Aside from identifying antigens for diagnostic purposes, proteomic studies have made invaluable contributions in drug target identification and understanding host-parasite relationships. This work has culminated in the identification of possible drug and vaccine targets for a range of parasites. In contrast, recent proteomic analyses on G. spinigerum aL3Gs have only focused on identifying immunoreactive proteins 10,12 . To date, immunoproteomics is the only proteomic analysis that has been used to study G. spinigerum. Other proteomic studies have been limited by the lack of a database for G. spinigerum protein identification. Therefore, protein identification has mostly relied on using nucleotide sequences from other nematode species. Recently, next-generation sequencing (NGS) was performed to provide a transcriptomic dataset from G. spinigerum aL3Gs 4 . This G. spinigerum database was applied in this study to facilitate proteomic and immunoproteomic identification of aL3G proteins and antigens. The findings of this study should improve current knowledge of the aL3G proteome and may help to elucidate the molecular functions and biological processes occurring in the parasite to ultimately determine novel drug targets. In addition, these results have identified immunogenic proteins which could be useful for the improvement of vaccines and diagnostic assays for gnathostomiasis.

Results
Protein profile of G. spinigerum third-stage larvae. In this research, protein and antigen profiles of aL3G were explored (Fig. 1). Proteins from aL3Gs were separated by 12% gel electrophoresis. A Coomassie blue-stained gel is shown in Fig. 2. A total of 14 gel pieces were cut and in-gel digestion and mass spectrometry analysis were performed. Proteins were identified by searching against the in-house G. spinigerum transcriptomic database 4 . With a 95% confidence interval cut-off, 687 proteins were identified (Supplementary Dataset 1). These G. spinigerum proteins were semi-quantified based on the exponentially modified protein abundance index (emPAI). The 20 most abundant aL3G proteins are shown in Table 1. The highly abundant proteins had an emPAI value that ranged from 0.99 to 289.02 and molecular weight (MW) that ranged from 7.5 to 250.9 kDa. Among the proteins identified, actin-2 expression increased by 43-fold compared with actin-5c, the second most highly expressed protein. Similarly, other structure-related proteins such as myophilin and actin-1 were also highly expressed in aL3Gs. In addition, several proteases such as zinc metalloprotease nas-14 and matrix metalloproteinase-like protein were also abundant in aL3Gs. Metabolic enzymes, including glyceraldehyde-3-phosphate dehydrogenase and nucleoside diphosphate kinase, were also highly expressed. Furthermore, enzymes responsible for protein folding, such as peptidyl-prolyl-cis-trans-isomerase 3, were also present at high levels in aL3Gs.
To gain more understanding into the biological processes occurring in aL3Gs, gene ontology (GO) was used to classify all identified G. spinigerum proteins (Supplementary Table 1). A total of 123 protein classes were found in aL3Gs. The top 20 protein classes are listed in Table 2. Proteins with unknown GO accounted for approximately 50% of the total identified G. spinigerum proteins. The major annotated protein classes observed in aL3Gs included embryo development ending in birth or egg hatching (GO:0009792), determination of adult lifespan (GO:0008340), and oviposition (GO:0018991). Other protein classes were involved in biological processes related  www.nature.com/scientificreports/  www.nature.com/scientificreports/ to energy, movement, morphological development, and enzymes for protein structure conformational changes. In addition to GO classification by biological process terms, proteins essential for parasite survival and host immunity evasion were also considered. Proteins relating to oxidation-reduction, proteinase-protease inhibitors, structure-movement, and energy were identified in the aL3G protein profile. The top 10 most abundant G. spinigerum proteins involved in these four categories are shown in Fig. 3. Glutathione S-transferase (emPAI of 0.37), peroxiredoxin (emPAI of 0.19), and nucleoredoxin (emPAI of 0.14) were the major antioxidant proteins expressed in aL3Gs. The most prevalent proteases were zinc metalloproteinase nas-14 (emPAI of 1.37), matrix metalloproteinase-like protein (emPAI of 1.13), and cathepsin B-like cysteine proteinase 6 (emPAI of 0.14), while the most abundant protease inhibitors were metalloproteinase inhibitor tag-225 (emPAI of 0.78), serpin I2 (emPAI 0.67), and serpin B4 (emPAI of 0.25). Meanwhile, actin-2 (emPAI of 289.02), actin-5C (emPAI of 6.7), and myoglobin (emPAI of 5.48) were the most highly expressed proteins in the structure-movement group. Additionally, glyceraldehyde-3-phosphate dehydrogenase (emPAI of 2.13), malate dehydrogenase, cytoplasmic (emPAI of 0.61), and phosphoglycerate kinase (emPAI of 0.37) were the main energy-related proteins expressed in aL3Gs.
Comparison of aL3G and human protein sequences. To identify potential G. spinigerum drug target candidates, the nucleotide sequences of the 10 most abundant proteins relating to oxidation-reduction, proteinase-protease inhibitors, and structure-movement were retrieved from the aL3G transcriptome. The blastx algorithm was used to search the lowest E-value for human proteins in the GenBank database compared with G. spinigerum sequences. The percent homology between G. spinigerum and human sequences is shown in Table 3. According to the alignment results, 10 G. spinigerum proteins including thioredoxin domain-containing protein 15, glutathione transferase omega-2, nucleoredoxin, glutathione S-transferase 1, carboxypeptidase A2, cathepsin B-like cysteinase 6, collagenase 3, myophilin, rho GDP-dissociation inhibitor 1, and imidazolonepropionase showed less than 50% homology to human sequences. Furthermore, no homology to human proteins was found in the following G. spinigerum proteins: chymotrypsin/elastase isoinhibitors 2 and 5; serpin B4; serpin I2; metalloproteinase inhibitor tag-225; zinc metalloproteinase nas-14; myoglobin; cuticle collagen 34; cuticle collagen 36; collagen alpha-1 (XVII) chain; and globin-like protein 9. Therefore, these G. spinigerum proteins might be potential candidates for drug development.
Protein identification in 24-kDa gel bands. Because identification of the 24-kDa aL3G crude worm antigen is an accepted diagnostic assay for gnathostomiasis, the aL3G proteins were separated by 12% gel electrophoresis then electro-transferred onto a membrane. Western blot analysis was performed using the sera of five individual patients with a confirmed diagnosis of gnathostomiasis as primary antibodies (Fig. 4). The aL3G proteins were transferred onto a membrane. The membrane was sliced into five 3 mm-strips. Each membrane strip (1)(2)(3)(4)(5) was individually incubated with sera from 5 different individuals diagnosed with gnathostomiasis As expected, the patient sera reacted with gel bands of approximately 24 kDa, a result reported to have high sensitivity and specificity for G. spinigerum diagnosis. Gel band numbers 9, 10, and 11 were excised for protein www.nature.com/scientificreports/ identification by mass spectrometry. The mass spectrometry data from three of the gel pieces were separately searched against NCBI and transcriptomic databases to identify their protein components (Table 4). Using the NCBI database, proteins from gel section numbers 9, 10, and 11 were identified as matrix metalloproteinaselike protein, unknown, and cyclophilin, respectively. Using the G. spinigerum transcriptomic database, matrix metalloproteinase-like protein and cyclophilin were identified from gels 9 and 11, respectively. Furthermore, 37, 26, and 23 proteins were also identified on gels 9, 10, and 11 using our database (Table 4). These data indicate the presence of other intriguing candidates, including antioxidative enzymes, cuticle collagens, and major antigens, which warrant further investigation as potential diagnostic targets.
2D-immunoblot aL3G antigen profile. To explore the aL3G immunome, aL3G proteins were separated by 2-DE and transferred onto a nitrocellulose membrane and later analyzed by western blotting (Fig. 5). Pooled patient sera served as the primary antibody. A total of 24 immunoreactive spots was observed. These spots were identified by comparing mass spectrometry data with the G. spinigerum transcriptomic database. A total of 115 proteins were identified as G. spinigerum antigens (Table 5). Spots 8, 17, and 22 demonstrated a lot number of protein identification at 21, 17 and 12, respectively. The proteins 32-kDa beta-galactoside-binding lectin, alkylated DNA repair protein alkB homolog 8, and mitogen-activated protein kinase were identified with the highest score in spot 8. Polyphosphoinositide phosphatase, UPF0378 protein, and heat shock factor binding protein 1 were observed with the highest confidence in gel 17. Matrix metalloproteinase-like protein, e3 ubiquitinprotein ligase trim13, and isocitrate isopropylmalate dehydrogenase domain containing protein were found in gel 22.
Potentially immunogenic proteins are known to be found in the 24-kDa region. A total of three, five, and 12 24-kDa proteins were identified in gels 20, 21, and 22, respectively. Interestingly, matrix metalloproteinase-like protein) was located in spots 20, 21, and 22, while cathepsin B was observed in spot 22 only. These two proteins are considered to be antigens in several parasites. Therefore, we predicted their protein properties, including molecular weight, isoelectric point, secretory and transmembrane domains, n-glycosylation, and o-glycosylation ( Table 6). The results showed that matrix metalloproteinase-like protein and cathepsin B were secreted proteins through signaling peptides and both of them contained n-glycosylation and o-glycosylation sites. As some clinical manifestations of G. spinigerum infection similar to Angiostrongylus spp., Strongyloides spp., and Sparganum spp. infection for example eosinophilic meningitis, larval migration and migration to the brain, distinguishing between G. spinigerum infection and the relating parasites might be helpful for diagnosis. Comparison of G. spinigerum matrix metalloproteinase-like protein and cathepsin B protein sequences with the highest % homology with protein sequences from humans and from related helminths, including Angiostrongylus spp., Strongyloides spp., and Sparganum spp., was performed. Matrix metalloproteinase-like protein and cathepsin B protein Figure 3. The top-10 most abundant proteins classified by molecular function relating to oxidation-reduction (redox), protease-protease inhibitor activity, structure-movenent, and energy metabolism. Their abundance was estimated semi-quantitatively by emPAI. In each category, the x-axis indicates the emPAI value and the y-axis specifies the protein identification. www.nature.com/scientificreports/ demonstrated 28%-31% and 34%-64% homology to related helminths and humans, respectively. Therefore, these proteins could be potential candidates for novel diagnostic assays to detect gnathostomiasis.

Discussion
Our data showed that actin-2 is highly expressed in aL3Gs. In addition, high levels of other G. spinigerum structure-related proteins, such as myophilin and actin-1, were also found. Actin is a cytoskeletal protein found in myofibrils in muscle cells and microfilaments in a variety of other cell types 13 . In trematodes, actin has been widely researched. It appears to be a crucial component of the platyhelminth tegument and is engaged in a variety of critical functions, including movements of secretory vesicles, muscle contraction, cytokinesis, and maintenance of cell shape 14 . In parasitic nematodes, cytoskeletal proteins play roles similar to those in trematodes and also have a function in nutrient uptake through transcuticular absorption 15 . The possibility of anthelmintic medications that attack helminth cytoskeletal proteins has piqued researchers' interest. In recent years, research on the microtubular network and its interaction with benzimidazole anthelmintics has been published 16,17 . Therefore, G. spinigerum structure-related proteins could be anthelmintic drug targets. In our study, GO classification of all identified aL3Gs proteins highlighted embryo development ending in birth or egg hatching (GO:0009792) as Table 3. List of Abundant proteins of aL3Gs classified by Gene Ontology (GO) into four groups under the molecular function category that are aligned with human proteins to identify unique proteins of aL3Gs. Chymotrypsin/elastase isoinhibitors 2 to 5* ---- Zinc metalloproteinase nas-14* ---- As well as GO classification by biological process terms, various protein classes required for parasite survival and host immunity evasion were also considered. The host immune system's oxidative burst can prevent intracellular parasite invasion and proliferation. Therefore, an antioxidant response promotes parasite survival, reduces inflammation, and changes the metabolism of the host cell 21 . In aL3Gs, the most abundant proteins relating to oxidation-reduction were glutathione S-transferase and peroxiredoxin. Moreover, these two proteins were also identified as G. spinigerum antigens (gels 9 and 10 in the 24-kDa region). Helminth glutathione S-transferase represents the main mechanism of detoxifying reactive oxygen intermediates produced by the parasite's endogenous metabolism or by the host immune system 22 . A 28-kDa glutathione-S-transferase was identified as an antischistosome vaccine candidate and recently evaluated in human clinical trials 23 . The glutathione-S-transferases have also been identified as potential immunotherapy or chemotherapy targets to treat infection with several other parasites, including Heligmosomoides polygyrus, Onchocerca gutturosa, and Dirofilaria immitis 24 . Peroxiredoxins are cysteine-dependent peroxidases that also play an important role in antioxidant, regulatory, and signaling systems. They are involved in defense against both endogenous and host-derived reactive oxygen species 25 . Treatment of mouse peritoneal macrophages with recombinant Toxoplasma gondii peroxiredoxin resulted in the production of IL-12p40 and IL-6 26 . This result suggests that peroxiredoxin induces both humoral and cellular immunological responses. Therefore, it could be used as a toxoplasmosis vaccine antigen. Peroxiredoxins have also been found in several parasitic helminths that affect humans such as Brugia malayi 27 , Onchocerca volvulus 28 , Taenia solium 29 , Opisthorchis viverrini 30 , and Schistosoma japonicum 31 . Furthermore, these proteins have been reported to be attractive therapeutic targets and vaccine candidates to treat and prevent helminth infections 32 . Nucleoredoxin was one of most abundantly expressed G. spinigerum proteins relating to oxidation-reduction. This protein is a newly discovered member of the thioredoxin family and plays a role in cell proliferation and differentiation 33 . Nucleoredoxin also regulates phosphofructokinase activity to balance between glycolysis and the pentose phosphate pathway 34 . However, there is limited information about this protein's role in parasites. Therefore, the function of nucleoredoxin in G. spinigerum needs to be further explored.

Structure-movement
Proteases and protease inhibitors are another important protein class expressed in G. spinigerum. The most prevalent proteases in aL3Gs were zinc metalloproteinase nas-14, matrix metalloproteinase-like protein, and cathepsin B-like cysteine proteinase 6. These three proteases had limited homology with human proteins. Therefore, they might be potential candidates for anthelmintic drug development. Moreover, they were also identified as G. spinigerum antigens (gels 9 and 11 in the 24-kDa region).
The ability of nematodes to molt is critical for their survival and development. Zinc metalloproteases, leucine aminopeptidases, and cysteine proteases are implicated in molting 35,36 . In the parasitic nematode Haemonchus contortus, the zinc metalloprotease nas-33 is required for molting and survival 37 . In Strongyloides ratti, collagenase plays a role in adult females at the time of migration through the host intestinal mucosa during oviposition 38    www.nature.com/scientificreports/ www.nature.com/scientificreports/ Owing to the various protease functions in parasitic worms, nematode proteases have been proposed as new anthelmintic and vaccine targets. The most abundant protease inhibitors in aL3Gs were metalloproteinase inhibitor tag-225, serpin I2, and serpin B4. There was no homology between these three protease inhibitors and human proteins. G. spinigerum serpins have been identified in excretory-secretory products of aL3Gs and have been shown to react with sera from gnathostomiasis patients 4 . Therefore, they may be good candidates for therapeutic and diagnostic development. Parasite protease inhibitors create a safer environment in the host by suppressing host protease activity and modulating host immunity 40 . Serpins from the parasitic worms Brugia malayi 41 , Ancylostoma caninum 42 , and Trichuris suis 43 have been identified. Brugia malayi serpin functions by neutralizing the immunostimulatory properties of the host cathepsin G 44 . Onchocerca volvulus serpin reduces the enzymatic activity of a panel of serine proteases including host elastase, chymotrypsin, trypsin, and cathepsin G 45 . Helminth protease inhibitors protect the worms from the hydrolytic action of host proteases and allow them to break through protective barriers and evade immunological responses. Protease inhibitors are, therefore, beneficial to parasite survival and colonization in the host 46 . Consequently, G. spinigerum proteases and protease inhibitors might have similar protective roles and may be good anthelmintic candidates.
Several G. spinigerum cuticle collagens-cuticle collagen 34, cuticle collagen 36, and collagen alpha-1 (XVII) chain-were identified in aL3Gs. These cuticle collagens had no homology with human proteins. However, these three cuticle collagens were conserved among other nematodes (data not shown). As a result, they may be suitable universal vaccine or drug targets for novel pan-nematicide therapeutics. Moreover, cuticle collagen 14 was identified as a 24-kDa G. spinigerum antigen and it showed less homology with protein sequences from other nematodes (data not shown). Therefore, this protein might be a potential target for the development of novel gnathostomiasis diagnostics and vaccines.
The exoskeleton of Caenorhabditis elegans, a free-living nematode, comprises a complex collagen matrix. Individual cuticle collagen gene mutations can result in exoskeletal abnormalities that change the shape of a nematode 47 . RNA interference targeted to the cuticle collagen genes of the root-knot nematode Meloidogyne incognita caused a 30.80%-35.00% reduction in the number of adult females and a 76.47%-82.59% reduction in the number of eggs. This study demonstrates the role of cuticle collagen genes in the structure and development of the nematode cuticle 48 . Furthermore, external enzymes that induce structural damage to the cuticle, such as papain, bromelain, collagenase, chitinase, and lipase, cause parasitic worms to die. After incubating Heligmosomoides polygyrus, Trichuris muris, and Protospirula muricola with plant cysteine proteases, the nematodes died as a result of cuticle damage 49 . These data suggest that cuticle collagens could be fascinating targets for nematicide development.
Western blot analysis has recently identified diagnostic targets for gnathostomiasis by detecting specific total IgG against a 24-kDa protein in aL3G extract 5 . This study detected specific IgG subclasses as well as total IgG against a 24-kDa antigen. When compared with other subclasses, IgG4 had the best sensitivity and specificity (91.6% and 87.8%, respectively) 6 . Although the detection of specific IgG against crude worm antigen (CWA) has a high sensitivity and specificity, the antigen preparation procedure is complex, time-consuming and laborious,  www.nature.com/scientificreports/ and batch-to-batch quality is inconsistent. Furthermore, natural sources of G. spinigerum are limited and depend on the season and climatic conditions. Therefore, identification of G. spinigerum candidate antigens is crucial for recombinant protein production to improve gnathostomiasis diagnosis. The 24-kDa diagnostic protein was first identified as matrix metalloprotease 9 . In another report, matrix metalloproteinase-like protein and cyclophilin with approximate molecular weights of 23-24 kDa were identified as G. spinigerum antigens 50 . However, identification of the protein composition of aL3Gs has been hampered by the lack of genomic sequence information required for proteomic analysis. In our study, the global aL3G antigen profile and specific 24-kDa aL3G antigens were identified using next-generation sequencing transcriptomic information 4 . Matrix metalloproteinase-like protein and cyclophilin were identified in the 24-kDa region, corresponding with previous reports. However, additional G. spinigerum antigen candidates were also identified. A total of 74 proteins were observed in 24-kDa bands on our SDS-PAGE gels. Our data suggest that glutathione S-transferase 1, cuticle collagen 14, major antigen, zinc metalloproteinase nas-4, major egg antigen, peroxiredoxin, and superoxide dismutase [Cu-Zn] might also be good candidates to explore for human gnathostomiasis diagnostic assays. On 2-DE, 13 proteins were also observed in 24-kDa bands.
A recombinant G. spinigerum matrix metalloproteinase protein has been produced. This protein was shown to have a sensitivity and specificity of 100% in a study comparing 40 patients with 30 healthy controls 51 . In another study, recombinant G. spinigerum matrix metalloproteinase protein was analyzed by immunoblot of serum samples from proven and clinically suspected cases of gnathostomiasis, patients with other parasitic diseases, and healthy volunteers. The sensitivity and specificity in this study were 100% and 94.7%, respectively 52 . Although matrix metalloproteinase protein is a potential candidate for G. spinigerum diagnostics, it may be important to  www.nature.com/scientificreports/ assess the different isoforms of this protein to improve its diagnostic capability. Our data also show that cathepsin B may be another good diagnostic candidate for human gnathostomiasis.
In this study, SDS-PAGE identified more proteins than 2-DE, possibly because 2-DE cannot resolve proteins that are too basic or too acidic, too large or too small. In addition, highly hydrophobic proteins are hard to dissolve in the 2-DE buffer system 53 . Therefore, integrating SDS-PAGE and 2-DE data from the 24-kDa region might be a useful way to explore novel diagnostic candidates for G. spinigerum. This study identified both matrix metalloproteinase like protein and cathepsin B from SDS-PAGE and 2-DE. Furthermore, these proteins were predicted to be secreted proteins and, therefore, might be useful for circulating antigen detection in patient serum. Moreover, both antigens were predicted to have glycosylation. Because glycosylation contributes to protein antigenic properties 54 , this post translational modification might need to be taken into account during any recombinant antigen production. the limitation of proteomics study using other G. spinigerum stages is to obtain adequate amount of parasite samples. In conclusion, the data from this proteomic and immunoproteomic analysis could help researchers better understand not only the parasite's biology but also possible targets for future treatments, diagnostics, and vaccines.

Methods
Preparation of third-stage G. spinigerum larvae. The aL3Gs were obtained from the livers of naturally-infected eels using an acid-pepsin digestion technique. Eel livers were brought from markets in Bangkok. They were chopped and subjected to 1% acid-pepsin at 37 °C for 2 h in a water bath stirred frequently, then were washed with tap water several times. Worms were then collected via a simple sedimentation technique. The collected worms were further identified using a dissecting microscope then washed again with tap water followed by normal saline solution (0.85% NaCl). Specie of G. spinigerum was confirmed by morphology. Approximately 20 collected worms were pooled in a microfuge tube and kept in -80 °C prior to analysis.

Sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE). The pooled aL3Gs
were ground using a mortar and pestle after short snap freezing. Lysis buffer was then added and further homogenization was performed using an ultrasonicator. After sonication, the mixture was centrifuged at 12 000 rpm for 15 min at 4 °C. The supernatant was collected and the protein concentration was determined by the Bradford method. Proteins from the lysate were separated by 12% SDS-PAGE gel and later stained with Coomassie Brilliant Blue G250 solution (Bio-Rad, Hercules, CA, USA). After running the SDS-PAGE, the gel was cut into 14 rectangles followed by in-gel digestion.
Two dimensional gel electrophoresis (2-DE). For the first dimension, a non-linear immobilized pH gradient (IPG) strip (pH 3-10; Amersham Bioscience, USA) was rehydrated overnight in immobilized pH gra- Immunoblotting. All methods were carried out in accordance with relevant guidelines and regulations.
All experimental protocols were approved by the Ethics Committee of the Faculty of Tropical Medicine, Mahidol University (MUTM 2020-058-02). Informed consent was obtained from all subjects and/or their legal guardian(s). Proteins from SDS-PAGE were transferred onto a nitrocellulose membrane. The membrane was cut into five strips and immunoblotting was performed using sera from five different patients as the primary antibody. The proteins separated by 2-DE were transferred onto a nitrocellulose membrane and pooled patient sera was used as the primary antibody. Membranes were blocked using 5% (w/v) non-fat milk in phosphate buffered saline (PBS) for 2 h at room temperature, then rinsed with PBS containing 0.05% (v/v) Tween-20. Serum samples diluted 1:200 in PBS containing 1% non-fat milk were added to the membranes and incubated overnight at 4 °C. After incubation, the membranes were washed three times with PBS containing 0.05% (v/v) Tween-200. Then horseradish peroxidase-conjugated goat anti-human IgG secondary antibodies were added and incubated for 1 h. Immunogen spots were visualized by detection of peroxidase activity using the Ultra TMB-Blotting Solution (ThermoFisher Scientific, UK). Immunoreactive protein spots were excised from silver-stained 2-DE gels and subjected to in-gel digestion.
In-gel tryptic digestion. Gel slices from SDS-PAGE and gels from 2-DE were destained until colorless. The former was destained with 50% acetonitrile (ACN, Sigma-Aldrich) in 50 mM ammonium bicarbonate (ABC, Sigma-Aldrich), while 30 mM potassium ferricyanide (K 3 Fe(CN) 6 , Sigma-Aldrich) and 100 nM sodium thiosulfate (Na 2 S 2 O 3 , Sigma-Aldrich) was used for the latter. After destaining, gel pieces were incubated in 4 mM DL-Dithiothreitol (DTT, Sigma-Aldrich) at 60 °C for 15 min. The DTT was subsequently removed and proteins were alkylated by adding 250 mM iodoacetamide (ICH 2 CONH 2 , Sigma-Aldrich) and incubated at room temperature in the dark for 30 min. The reaction was quenched with 4 mM DTT and dehydrated in 100% ACN. To digest the proteins, the gel pieces were again rehydrated with 10 ng/µL trypsin in 50 mM ABC and incubated at 37 °C overnight. The peptides were recovered by adding ACN. The supernatant was collected and dried using a vacuum centrifuge (TOMY, Japan). Dried peptides were resuspended in 0.1% formic acid for LC-MS/MS analysis.

Mass spectrometry.
A MicroTOF Q II mass spectrometer interfaced with an Ultimate™ 3000 nano-LC system was used for LC-MS/MS. An Acclaim PepMap RSLC 75 µm × 15 cm nanoviper C18 column with 2 µm particle size and 100 Å pore size (Thermo Scientific, Waltham, MA) was used. Data were acquired using a MicroTOF Q II mass spectrometer set with a scan range of 500-3500 m/z. Mass spectrometry data were analyzed using the MASCOT search engine 2.3 (Matrix Science, Ltd) for peptide identification and a previously reported in-house transcriptomic database 1 as the reference proteome. Search parameters were set as follows: one miss cleavage; trypsin digestion; 0.8 Da peptide tolerance; ± 0.8 fragment mass tolerance; carbamidomethyl (C) and oxidation (M) variable modifications. Significance threshold was 0.05. The protein abundance was also determined based on the LC-MS/MS output 55 .
Bioinformatic analysis. After protein identification using the in-house transcriptomic database based on the recently published ES proteome of infective stage G. spinigerum larvae 4 , the identified proteins were further classified using the GO database (http:// www. geneo ntolo gy. org) to determine and predict the biological processes affected by these parasitic proteins. Furthermore, the sequences of proteins identified by the transcriptomic data were also subjected to the Basic Local Alignment Search Tool (BLAST) translated (BLAST:blastx) and searched against the human non-redundant protein database to identify possible drug and vaccine targets. The % identity, E-value, and query coverage were reported. For matrix metalloproteinase-like protein and cathepsin B, property predictions and molecular weight and isoelectric point calculations were performed by the Expasy Compute pI/Mw tool (https:// web. expasy. org/ compu te_ pi/). Prediction of signaling peptides, non-classical secretory domains, transmembrane protein domains, n-glycosylation, and o-glycosylation were performed using SignalP version 5.1 (https:// servi ces. healt htech. dtu. dk/ servi ce. php? Signa lP-5.0), SecretomeP version 2.0 Server (https:// servi ces. healt htech. dtu. dk/ servi ce. php? Secre tomeP-2.0), TMHMM version 2.0 (https:// servi ces. healt htech. dtu. dk/ servi ce. php? TMHMM-2.0), NetNGlyc version 1.0 (https:// servi ces. healt htech. dtu. dk/ servi ce. php? NetNG lyc-1.0), and NetOGlyc version 4.0 (https:// servi ces. healt htech. dtu. dk/ servi ce. php? NetOG lyc-4.0), respectively. For sequence alignment, all sequences were retrieved from the non-redundant protein sequence NCBI database. The alignments and identity calculations were performed using the Clustal Omega software.

Dat availability
The dataset used in this study might be shared upon reasonable request to Onrapak Reamtong, PhD.